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Abstract 

In this paper we propose a strategy to approximate incompressible hydrostatic 
free snrface Enler and Navier-Stokes models. The main advantage of the proposed 
models is that the water depth is a dynamical variable of the system and hence 
the model is formnlated over a fixed domain. 

The proposed strategy extends previons works approximating the Enler and 
Navier-Stokes systems nsing a mnltilayer description. Here, the needed closnre 
relations are obtained nsing an energy-based optimality criterion instead of an 
asymptotic expansion. Moreover, the layer-averaged description is snccessfnlly 
applied to the Navier-Stokes system with a general form of the Canchy stress 
tensor. 

Keywords : Incompressible Navier-Stokes equations, incompressible Euler equations, 

free surface flows, newtonian fluids, complex rheology 
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1 Introduction 

Due to computational issues associated with the free surface Navier-Stokes or Euler 
equations, the simulations of geophysical flows are often carried out with shallow water 
type models of reduced complexity. Indeed, for vertically averaged models such as the 
Saint-Venant system j7|, efficient and robust numerical techniques (relaxation schemes 
|9], kinetic schemes |25l E],. ..) are available and avoid to deal with moving meshes. 

In order to describe and simulate complex flows where the velocity field cannot be 
approximated by its vertical mean, multilayer models have been developed PE] SI El 
mm Unfortunately these models are physically relevant for non miscible fluids. 

In [TBl El El ES], some authors have proposed a simpler and more general formulation 
for multilayer model with mass exchanges between the layers. The obtained model 
has the form of a conservation law with source terms, its hyperbolicity remains an open 
question. Notice that in j5] the hydrostatic Navier-Stokes equations with variable density 
is tackled and in |26) the approximation of the non-hydrostatic terms in the multilayer 
context is studied. With respect to commonly used Navier-Stokes solvers, the appealing 
features of the proposed multilayer approach are the easy handling of the free surface, 
which does not require moving meshes (e.g. ini). and the possibility to take advantage 
of robust and accurate numerical techniques developed in extensive amount for classical 
one-layer Saint Venant equations. Recently, the multilayer model developed in [16] has 
been adapted in |15] in the case of the /i(I)-rheology through an asymptotic analysis. 

The objective of the paper is twofold. First we want to present another derivation 
of the models proposed in [HI E] [26], no more based on an asymptotic expansion but 
on an energy-based optimality criterion. Such a strategy is widely used in the kinetic 
framework to obtain kinetic descriptions e.g. of conservations laws [201125] . Second, we 
intend to obtain a multilayer formulation of the Navier-Stokes system with a rheology 
more complex than the one arising when considering newtonian fluids. 
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The paper is organized as follows. In Section we recall the incompressible hydro¬ 
static Navier-Stokes eqnations with free surface with the associated boundary conditions. 
In Section]^ we detail the layer averaging process for the Euler system and obtained the 
required closure relations. The proposed layer-averaged Euler system is given in Sec¬ 
tion]^ and its extension to the Navier-Stokes system with a general rheology is presented 
in Section m 


2 The Navier-Stokes system 

We consider the two-dimensional hydrostatic Navier-Stokes system izu describing a free 
surface gravitational flow moving over a bottom topography Zb{x). For free surface flows, 
the hydrostatic assumption consists in neglecting the vertical acceleration, see [TUI [TUI lUU] 
for justifications of the obtained models. 


2.1 The hydrostatic Navier-Stokes system 

We denote with x and the horizontal and vertical directions, respectively. The system 
has the form; 

( 1 ) 
( 2 ) 
(3) 

and we consider solutions of the equations for 


ou aw 
dx~^ dz 

du dv? duw dp dT^^x 

dt ^ dx ^ dz dx dx 

dp 

dz 


+ 


as. 


dz ’ 


= -9 + 


as. 


dx 


+ 


as. 


dz ’ 


t > to, X G M, Zb{x) < z < p{x,t), 

where p{x, t) represents the free surface elevation, u = {u, w)'^ the velocity vector, p the 
fluid pressure and g the gravity acceleration. The water depth is H = p — Zb, see Fig. 
The Cauchy stress tensor S-p is defined by S^^ = —pid + S with 



and S represents the fluid rheology. 

As in Ref. ini, we introduce the indicator function for the fluid region 


^{x,z,t) 


1 for (x, z) & Vt = {{x, z)\zb < z < p}, 
0 otherwise. 


(4) 


The fluid region is advected by the flow, which can be expressed, thanks to the incom¬ 
pressibility condition, by the relation 


dip difu dipw 
dt dx ^ dz 


0 . 


(5) 
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Figure 1; Flow domain with water height H{x,t), free surface ri{x,t) and bottom Zb{x). 


The solution ip of this equation takes the values 0 and 1 only but it needs not be of 
the form (|^ at all times. The analysis below is limited to the conditions where this 
form is preserved. For a more complete presentation of the Navier-Stokes system and 
its closure, the reader can refer to |21| . 

Remark 2.1 Notice that in the fluid domain, Eq. I® reduces to the divergence free 
condition whereas across the upper and lower boundaries it gives the kinematic boundary 
conditions defined in the following. 


2.2 Boundary conditions 

The system 0-0 is completed with boundary conditions. We not consider here lateral 
boundary conditions that can be usual usual inflow and outflow boundary conditions. 
The outward unit normal vector to the free surface and the upward unit normal 
vector to the bottom are given by 


n, = 


i+(i) 


drj 

dx 


Hfo = 


1 + 


\dx ) 


dzb 

dx 

1 


-Sb 

Cb 


respectively. We use here the same definition for Sb{x) and Cb{x) as in [9], Cb(x) > 0 is 
the cosine of the angle between and the vertical. 


2.2.1 Free surface conditions 

At the free surface we have the kinematic boundary condition 

drj dp 

— +Us^ -Ws = 0, 

ot ox 


( 6 ) 


where the subscript s indicates the value of the considered quantity at the free surface. 
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Assuming negligible the air viscosity, the continuity of stresses at the free boundary 
imposes 

(7) 

where p°‘ = p°'{x,t) is a given function corresponding to the atmospheric pressure. 
Within this paper, we consider = 0. 


2.2.2 Bottom conditions 


The kinematic boundary condition at the bottom consists in a classical no-penetration 
condition; 

n n 

Ufc ■ nf, = 0, or Ub- -tCfe = 0. 

ox 

For the stresses at the bottom we consider a wall law under the form 


Ernfo - (nft ■ T,Tnb)nb = KUb (9) 

and for t^ =*(cb, Sb), using (|^ we have 

tb ■ Sj-nfe = —Ub, (10) 

Cb 

If K(ub, H) is constant then we recover a Navier friction condition as in [1^. Introducing 
a laminar friction ki and a turbulent friction kt, we use the expression 


K(ub,i 7 ) = ki + ktH\uk\, 

corresponding to the boundary condition used in Another form of fi:(ub, H) is used 
in [9], and for other wall laws the reader can also refer to Due to thermo-mechanical 
considerations, in the sequel we will suppose ^(ub, H) > 0, and ^(ub, H) will be often 
simply denoted by k. 


2.3 Other writing 

For reasons that will appear later, we rewrite 0-0 under the form 


du dw 
dx dz 

du du^ duw dr] dTj^x dT^^z d'^ 

dt dx dz ^ dx dx dz dx‘^ 


'^zxdz\ 


d^zz 

dx 


( 11 ) 

( 12 ) 


where Eq. (12) has been obtained as follows. Integrating Eq. (|^ from z to p and taking 
into account the boundary condition ([^ gives 


d n 

p = g{p - z) - — / T.^xdzi S, 


( 13 ) 


Inserting the previous expression for p in Eq. ([^ gives Eq. ( [I^ . 
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2.4 Energy balance 

Lemma 2.2 We recall the fundamental stability property related to the fact that the 
hydrostatic Navier-Stokes system admits an energy that can he written under the form 


d 


d 


cr] 


— E dz + — 


dt 


with 




dx 


'Zb 


A 

dx 




'Zb 


\dx 


dz 


dx 


E = ^ 
2 


gz. 


1 j 

dz 

n ^ 2 

dz - 

(14) 
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(15) 


Proof of lemma 12.21 The way the energy balance (14) is obtained is classical. Con¬ 
sidering smooth solutions, first we multiply Eg. ^ by u and Eg. (|^ w then we sum 
the two obtained eguations. After simple manipulations and using the kinematic and 
dynamic boundary conditions (g-(§ , we obtain the relation 


d 


d 


— E dz — 


dt 


'Zb 


dx 


'Zb 


_ 

/ ^xx o dz 

Jzb dx 


u{E + p) - uT^xx - dz 

du , n dw . 


TjT.,—dz — / -—Yi,T.dz — 


'Zb 


'dz 


'Zb 


dx 


'Zb 


\ ^dz-Zu^ 


By using Eg. 0 and replacing p by its expression given by ( [I^ in the previous relation 
gives the result. 


3 Depth-averaged solutions of the Euler system 


In this section, neglecting the viscous effects in Eqs. 0-0 , we consider the free surface 
hydrostatic Euler equations written in a conservative form 


d(p difu difw 

dt dx dz 

dtpu dg)u^ dtpuw dp 

dt ^ dx ^ dz ^ dx 


dp 

dz 


= -T9, 


(16) 

(17) 

(18) 


with g) deffned by (|^. 
and ([^ that reduces to 


This system is completed with the boundary conditions ([^,(|^ 


Ps = 0 . 


(19) 
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From Eqs. (18),(19), we get 

p = ipg{Ti-z). (20) 

The energy balance associated with the hydrostatic Euler system is given by 


d 


d 


I I '^{E + P)dz = 0, 


Zb 


dx 


V 

Zb 


( 21 ) 


with E dehned by (15). 


3.1 Vertical discretization of the fluid domain 

The interval [zb,ri] is divided into N layers {La}ae{i,...,N} of thickness laH{x,t) where 
each layer La corresponds to the points satisfying z e La{x,t) =] 2 :a-i/ 2 , 2 ^a+i/ 2 [ with 


r Za+i/2{x,t) = Zb{x) + 

\ ha{x, t) = Za+l/2{x, t) - Za-ll2{x, t) = laH{x, t), N}, 


with Ij > 0, 

We also dehne 


1, see Fig. 


^a+l/2 + Za-1/2 ha r , 

^- = Za-1/2 + Y^ a = {1,. . . , A/}. 

We hnally introduced the distance between the midpoints of the layers, 

ha+i + ha 


ha+l/2 Za+1 Za 


a = {1,... ,iV - 1}. 


( 22 ) 


(23) 


(24) 



Figure 2: Notations for the multilayer approach. 
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3.2 Layer-averaging of the Euler solution 


In this section we take the vertical average of the Euler system and study the necessary 
closure relations for this system. 

Let us denote {f)a the integral along the vertical axis in the layer a of the quantity 
/ = fiz) he. 

{f)a{x,t)= / (25) 

Jr 

where lz&La(x,t){z) is the characteristic function of the layer a. 

The goal is to propose a new derivation of the so-called multilayer model with mass 
exchanges |6l [5] using the entropy-based moment closures proposed by Levermore in [19] 
for kinetic equations. This method has already been successfully used by some of the 
authors in m 

Taking into account the kinematic boundary conditions (|^ and ([^, the layer-averaged 
form of the Euler system (16)-(18) writes 

d 


d 

~ Ga+1/2 — Ga-1/2, 
d d dp 

— {(pu)a + )“ “h ~ '^a+l/2Ga+l/2 ~ WQ-l/2Ga-l/25 

d d 

— -t- —{ipzu)a = {(pw)a + Za+ll2Ga+ll2 “ ^a-l/2Ga-l/2, 


(26) 

(27) 

(28) 
(29) 


for a G {1,..., N} and where p is defined by Eq. (20). The quantity Ga+ 1/2 is defined 
by 

^^a+1/2 , dZa+l /2 

+ Ua+l/2- 


Gr 


:+l /2 — 7 ’a-|-l /2 


dt 


dx 


w. 


0+1/2 I , 


(30) 


and corresponds to the mass flux leaving/entering the layer a through the interface 
Za+i/ 2 - The value of (^ 0 + 1/2 is equal to 1 for every a. Notice that the kinematic boundary 
conditions ([^ and (|^ can be written 

Gi/2 = 0, GAr+ 1/2 = 0. (31) 

These equations just express that there is no loss/supply of mass through the bottom 
and the free surface. Taking into account the condition (31), the sum for j = 1,... a of 

d 


the relations (26) gives 




d 


The quantities 


0+1/2 


u. 


dt 




/=! 


i=i 


(+ 1/2 — ^a+J/ 2 yi)y 


(32) 


(33) 
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corresponding to the velocities values on the interfaces will be dehned later. Notice that 
when using the expression (32), the velocities Wq,+i /2 no more appear in Eqs. (26)-(29) 
and thus need not be dehned. 


Equation (29) is a rewriting of 


dip dpu dpw 
dt dx ^ dz 


- 1/2 


, \ ! , dpw\ 


using again the kinematic boundary conditions. Notice also that because of the hydro¬ 
static assumption, Eq. (29) is not a kinematic constraint over the velocity held but the 
dehnition of the vertical velocity {pw)a- The form of Eq. (29) is useful to derive energy 


balances but other equivalent writings can be used, see paragraph 4.2 


Simple manipulations allow to obtain the system (26)-(30) from the Euler sys¬ 
tem (16)-(18) with ([^ and (|^ e.g. for Eq. (26), starting from (16) we write 

^ d^ d^ _ 

^dt dx dz 


and using the Leibniz rule to permute the derivative and the integral directly gives (26). 
Likewise, the Leibniz rule written for the pressure p gives 


A = 

^dx'" 


.+1/2 dp d dz^+i/2 dZa-l/2 

— dz = -^{p}a -Pa+112 -^- Pa-1/2- 

' y~4i-Y> ' 


Za-l/2 


dx 


dx 


dx 


and from (28),(19), we get 


N 


Pa+l/2=P{x,Zc,+ l/2,t) = ^ {pg)j. 

j=a+l 


(34) 


From Eq. (20), we also have 

l‘^a + \/2 Q 


d-i) = 

^dx'" 




dzh 


^a-l/2 


Relation (20) also leads to 

P = Pa+l /2 + 9 P{Za+l /2 - Z) = Pa-l /2 + gp{Za-l /2 “ z), 


and hence 


/ \ / \ Pa+l/2 + Pa-l/2 , > , 9 / s2 

\P)a = {P)a -^- = {p)aPa+l/2 + a' 


Therefore, the system (26)-(30) can be rewritten under the form 
d d 

-^{p)a + “ Ga+l/2 — Gq-I/2, 

d d 

-^{pu)a + ^ + {p)a) = Ua+l/2Ga+l/2 — Ua-l/2Ga-l/2 


dx 


(35) 


( 36 ) 
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dz, 


Pa+l/2- 


'a+1/2 


dz„- 


dx 


Pa-l/2- 


01 - 1/2 


dx 


(37) 


d d 

— {(pz)a + -^{y:>Zu)a = {pw)a + 2;o+l/2G^o+l/2 ~ Za-1/2Ga-1/2, (38) 


with (34),(35) and completed with relations (32). 


Considering smooth solutions, multiplying (17) by u and integrating it over the layer 


a gives, after simple manipulations, the energy balance 

,2 


d 


d 


-{E)^ + -(n(E + p))^ = 


u 


0 + 1/2 


+ Pq+1/2 + gZa+l/2 1 G'q,+i/2 


^0-1/2 \ _ dZa + l /2 dZa - ll 2 , . 

+ Po-1/2 + 5 '^o- 1/2 ) Cq,_i /2 — Po+ 1/2 -- 1 “ Po-1/2--5 (39) 


dt 


where E = E{z]u) is defined by (15). The sum for a = 1,... ,iV of the relations (39) 
gives 

^ N ^ N 

^ + p))a = 0- 


dt 


a=l 


a=l 


Therefore the system 


completed (^32h, (34) and (35) has three equations 


with three unknowns, namely {^p)ai {pu)a and {ipw)a and closure relations are needed 
to define {ipu^)a, {pzu)a and u{x, Za+l/2,t). 


3.3 Closure relations 

If u'^ is defined as the deviation of u with respect to its layer-average over the layer a, 
then it comes for z & La 

PU = — -h LfU^, 


{l)c 

with {pu'^) = 0. Following the moment closure proposed by Levermore 
the minimization problem 

mm{{ipE{z;u)})a- 

< 

The energy E{z-,u) being quadratic with respect to u we notice that 

2 \ _ {Lu)1 , 2{ipuu')a , / / , 32 \ 


{ipU^)a = 


> 


{P)a 

Ml 

(l) 

(ipu) 


+ 


{l)c 




+ {pWaY) 


a. 

2 

a 


(T’). 


(40) 
we study 

(41) 


(42) 


Equation (42) means that the solution of the minimization problem (41) is given by 


{ipE z] 


{pu)c 

{P)a 


)a = m:in{{ipE{z]u)})c 


(43) 
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and 


{^E 


/ M 

V ’ {^)c 




Since the only choice leading to an equality in relation (42) corresponds to 

(ifu). 


ipu = 


{V)c 


for z e La, 


this allows to precise the closure relation associated to a minimal energy, namely 

{Ln)l 


{(fU )a = 


{L)a ’ 


{ipzu)a = {(pz), 


{lu)c 
^ (V?)a 


(44) 

(45) 

(46) 

(47) 


It remains to dehne the quantities Ua+i/ 2 - We adopt the dehnition 




0 + 1/2 


{ipu)o 
- } .(V’>c 


{tfu)a + l 

{Da + l 


if G, 
if G, 


0 + 1/2 ^ 


< 0 


0 + 1/2 


> 0 


(48) 


corresponding to an upwind dehnition, depending on the mass exchange sign between 
the layers a and a + 1. This choice is justihed by the form of energy balance in the 
following proposition. 

Proposition 3.1 The solutions of the Euler system (16)-(18) with (@,(0) satisfying the 
closure relations (46)-(48) are also solutions of the system 


d d 

~^{T)a + = Ga+1/2 ~ 1 ^ 0 - 1/25 

d , \ Of {(pu)l^ / \ \ ^ ^ 

— {ipU)a + ^ i ~(^)-^ ^ ~ '^ 0 + 1 / 20 ^ 0 + 1/2 ~ Ua-l/2Ga-l/2 


(49) 


dz. 


Pa+l/2- 


0 + 1/2 

dx 


d 


d 


777 


dt 


dx 


^ {T)a 


^^0-1/2 

'Pa-112 rs , 

dx 

(50) 

0+1/2 ~ ^0-1/2^0-1/25 

(51) 


completed with relation (32). The quantities {p)a andpa+ 1/2 ore defined by (34) and (35). 

This system is a layer-averaged approximation of the Euler system and admits - for 
smooth solutions - an energy equality under the form 


m 


N 


a=l 




(t) 


d y^,{ipu)c 


dx^^ (if). 




{t)c 


N 


1 y-v f {(pu)a+l 


Q:=l 


(t) 


a+1 


{tu)c 


1^0+1/21. (52) 
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Remark 3.2 Instead of (48), the definition 


u. 


a+l/2 


fip)a+l ( 9 ^) 0 : 


( 53 ) 


is also possible and gives a vanishing right hand side in (52). But such a choice does 


not allow to obtain an energy balance in the variable density case and does not give a 
maximum principle, at the discrete level, see m- Simple calculations show that any other 


choice that (|48|) or (|53|) leads to a non negative r.h.s. in (|52|), see Eg. (|54|) in the proof 
of prop. \3.1 


Remark 3.3 It is important to notice that whereas the solution H,u,w,p of the Euler 
system (16)-(19 ),(|^,@^ a/so satisfies the system (36)-(38), only the solutions H,u,w,p 
of the Euler system (16)-(19),([^, ( [^ s atisfying the closure relations ([4^-([4^,([4^ are 
also solutions of the system (49)- (|52[ ). On the contrary, any solutions {‘pu)a, 

{(pw)a and {p)a of (49)-(51) with (48) are also solutions of (36)-(39). 


Proof of prop. 3.1| Only the manipulations allowing to obtain (52) have to be detailed. 
Eor that purpose, we multiply (ISOl) by giving 

1——I y-P/oL 


/ d {(pu)a 
\dt {ip) a dx 



{y>u)c 


= U. 


a+l/2Ga+l/2 ~ ''^o-l/2<-Ja-l/2 




+ 


dz, 


'a+l/2 

dx 


dZry- 


Pa+1/2 


0 - 1/2 

dx 


Pq-1/2 


(/3^^)o 

(V^)« ’ 


and we rewrite each of the obtained terms. 

Considering first the left hand side of the preceding eguation excluding the pressure 
terms, we denote 


/ d {ipu)a d f {ipu)a 

\dt {ip)a dx V {p)a )) {p>)a 


and using (26) we have 


‘■u^a 


^ ( {^U)l \ df {ipu)c, {pu)l , 

dt \2{ip)a) dx V (V?)o 2{ip)a) 2{ip)l 



Now we consider the contribution of the pressure terms over the energy balance i. e. 


[ d{p)a dZa+112 , dZa-l/2\ {<pu) a 

\ dx ^“+ 1/2 Pq- 1/2 J 


Using (35) we get the eguality 


dza+il2 

Po-i- 1 / 2 —^- 

dx 


Pq-1/2 


dZa-ll2 

dx 


{p)a dhg 
{ip)a dx 


{9P)a 


dZg 

dx ’ 
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holds, it comes 


4.<. ({P>< 


is I 

I I <">' 


i^u), 

{^)a 

(ipu). 


ip).l 


{^)c 

(ipu) 


dx \ (if), 
{p)a d{ipu)a , d 


{p)a {^U)a d{p)a , , . Oz^ 

T ~\—-— +9{^u}c - 

ma ma OX 


dx 


+ TT aKz, 


{^)c 


{p)a dx dx 

{p)a d{pu)^ , ^ I ..u ^ 

/ \ r-v O 1 y'^Ci^C 

{ip}a OX ax 


i^u), 

' {^)a 

(pu). 


d 


- Za-^ {9{v>u)a) 

(^>w "4 «»*’>"> 


— 9Za {Ga+1/2 — Ga-1/2) 


= -iziip) 


dx 


{P'^)c 

{P)a 


{p)a d{pu)o, , d 


{ip) a dx dx 


+ 7^ 9haZ, 


(pu) 


d 


/ \ I + 7u i9haZa) 

{p)a } dt 


dz, 


~9ha^Z, - 9Zq {Ga+l/2 ~ Ga-l/ 2 ) 


dt 

Let us rewrite Ip^a under the form 

{pu)c 


4« ((4. 


{P)c 


d 


d 


+ ZZ idhaZa) + 9 ZT { Kz, 


dt 


dx 


{u)c 

hr. 


~9 4a+l/2GQ,+i/2 — Za-il2Ga-ll2) + Jp,cn 


with 


{p)a dx 


{p)r,d{pu)a r,h M (r MT \ 

~ ^” 2 ” \ya+l/2 + ^a-l/2) 


Since we have 


{p)a d{pu)a _ 

{p)a dx {p) 


I ^ 

isr^4_-\ /9 — CjTrv- 


we obtain 


j _ dza+ 1/2 

'Jp^a Pa-\-l/2 W, Pa.—ll2 


dt 

Then summing Iu,a und Ip^a gives 


I 


{Tu)r 

{T)a 


)« + Tizi 


d ,{pu)r 


dx {p)r 


E z; 


a+1/2 

dZa-ll2 

dt 

{pu) 


a- 1/2 


dha 

dt 


{t)c 


~ Pa+l/2Ga+l/2 + Pa-l/2 Gol-1/2- 


T {p)a ) )a 


U 


0 + 1 / 2 - 


{Pu)a _ 1 / {pu) 

{p)a 2 V {P)c 


G. 


0+1/2 


Ury — 


a-1/2- 


{Pu)a _1 / {pu) 
{p)c 2 V {P)c 


Gr,- 


0 - 1 / 2 - 


Finally, the sum of the preceding relations for a = 1,..., N 


® ^ f ^ I + (p>« I)«= 


dt 


a=l 

N 

E 

a=l 


dx^^ {p)c 


{t)c 


u 


0 + 1/2 


{pu)a _ {pu)a+i \ _1 / {pu) 

{P)a {P)a+1 ) 2 V {P)c 


^ 1 / (</?m)„+i 


{t) 


a+l 


G^o+1/25 (54) 


and the definition (48) gives relation (52) that completes the proof. Notice that any other 


choice than (48) or (53) leads to a non negative r.h.s. in (54), see remark 3.2 
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4 The proposed layer-averaged Euler system 


4.1 Formulation 


The closure relations (46)-(47) motivate the definition of piecewise constant approxima¬ 
tion of the variables u and w. 

Let us consider the space of piecewise constant functions defined by 

a e 

Using this formalism, the projection of u and w on P^)^ is a piecewise constant function 
defined by 


N 


X^{x,z,{za],t) = 5 ^ 1 ; 


^a-l/2)2a, + l/2 




(55) 


a=l 


for X G {u,w). In the following, we no more handle variables corresponding to vertical 
means of the solution of the Euler equations (16)-(18) and we adopt notations inherited 
from (55). 


By analogy with (49)-(51) we consider the following model 

dha d{haUa) _ 

dx ~ ’ 

a=l Q=1 

dh^u„ d , 2 , u \ ^ ^ 

+ ^ [n-aUa + iT'CtPa) — Ua+l/2^a+l/2 ~ Ua-ll2^a-ll2 


dt dx 


dz. 


a+l/2 


dz^_ 


d ( ^0+1/2 ^a-1/2 \ ^ I '^a+1/2 ^a-l/2 \ , 

7W -^- +7^ -^- Ua = haW, 


— Z 


d zi 


dx 
— zti 


Pa+l/2 


a-1/2 

dx 


Pa- 1 / 2 , 


dt 


dx 


by analogy with (32) 


d d 

Ga+l/2 = ’ 

i=i i=i 


and we have pa, Pa+ 1/2 given by 


(56) 


(57) 


+^a+l/2G'a+l/2 ~ ^0-1/21^0-1/2, (58) 


(59) 


hr 


N 


N 


Pa= 9 


Y 

j=a+l 


and Po+i/2 - g Y 

j=a+l 


( 60 ) 


The definition of Ua+ 1/2 is equivalent to (48) i.e. 


u 


0+1/2 


Ua if Go+1/2 < 0 

Mo +1 if ^0+1/2 > 0 
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The smooth solutions of (56)-(58) satisfy the energy balance 

0 0 f \ 

^ {Ua {Ea + haPa)) = ( 'WQ+l/2Wa-^ + Pa+1/2 + 5'^a+l/2 1 ^0+1/2 

“a \ 

i-l/2'Wa- — + Pa-ll2 + 5'^o-l/2 I bra-1/2 




■ Pa+1/2- 


a+1/2 

at 


+ Pa-1/2- 




a-1/2 


at 


( 61 ) 


with 




-^a — "1 “ + §(^a+l/2 ~ ^a- 1 / 2 ) — ^a ( ^ 


Adding the preceding relations for a = 1,..., A^, we obtain the global equality 


Ft 


N 


0 


N 


N 


^ ^ Eg I T I ^ ^ (T/q, T hoiPa') I ^ ^ ^ (^a+1/2 ^a) |f?a+l/2|' (62) 


ya=l 


ya = l 


a=l 


Using (60), the pressure terms in (57) can be rewritten under the form 

OZa-l/2 


0 OZa+1/2 

yhotPa) a Pa+1/2 T a 

Ox Ox Ox 


4.2 The vertical velocity 




Ozh 


The equation (58) is a definition of the vertical velocity given by (55). The quantities 
Wa are not unknowns of the problem but only output variables. Indeed, once El and 
have been calculated solving (56),(57) with (59), the vertical velocities Wa can be 
determined using (58). 


Using simple manipulations, Eq. (58) can be rewritten under several forms. In par¬ 
ticular, the following proposition holds 

Proposition 4.1 Let us introduce w = w{x,z,t) defined by 


Ou^ Ow 

^ = 0, 


(64) 


Ox Oz 

The quantity w is affine in z and discontinuous at each interface Za+\i 2 , w can he written: 

OUa 


w = kn, — z 


Ox ’ 


( 65 ) 


with ka = ka{x,t) recursively defined by 

0{zbUi) 


ki = 


Ox 


0 


ka+l kot T (.2a+l/2(^a+l ^a))- 
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Therefore we have 


f^a + l/2 


Wdz = haWa, 


( 66 ) 


' ^a-1/2 


meaning the quantities w is a natural and consistent affine extension of the layer-averaged 

e la'i 

duo 


quantities Wa defined by (58). Using (66), an integration along the layer a of (65) gives 

.v2 


or 


h W -h k ^a-l/2 dUg _ _ 

2 ox 


, dUa .. . 

Wa = ka- = w{Za) 

OX 


dx 


(67) 

( 68 ) 


Proof of prop. 4.1 A simple integration along z of equation (64) using (^8| gives 


9 r N, 

w = —— f u dz, 


(69) 


and therefore, for z E Li we get 


d 


d 


w 


I '>^idz = - — {{z-z,)u,), 


'Zb 


dx 


i.e. 


d . , dui 

W = -^{ZbUi) - z——. 

dx dx 


For z G La, relation (69) gives 

a—l 


W = 


d d 

H ^a-l/2)Wa), 


i=i 


and we easily obtain 


, dua 
w = ka- z——. 

dx 


Now we intend to prove (66). 


Using the definition (23), relation (58) also writes 
d , , ^ A d{hjUj) 


Q.— 1 


1 ^ / 1 \ uutjaj) 

haWa \ZahaUaj Za+1/2 / ^ fpr 1/2 / ^ 


dx 


i=i i=i 

leading to a new expression governing Wa under the form 


djhjUj) 

dx 


a—l 


ha d{haUa) , 


2 dx 


i=i 


dx 


+ haUa 


dZa 

dx 


(70) 


(71) 
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And from (70), we get 


ol—1 


r^O( + l/2 -^ ^ ^ dUfy 

/ wdz = -ha^-^{hjUj) + ha — {Za-l/ 2 Ua)-haZa 

^a-ll2 j = l 

a—1 


dx 


J~l p) py 


i=i 


corresponding to (71) and proving the result. 


Using also (63), we are able to rewrite the system (56)-(58) under the form 

N 


N 


E dha d{haUay _ 


1 1 
a=l a=l 

dhnUry d 


dx 


= -- 


2 dx 


a—1 

i=i 


dx 


+ Uc 


^ {haU^a + ~poiH^ — —gha-^ + Ua+l/2Ga+l/2 ~ Wa-l/2Ga-l/2, (73) 

1 d{haUa) ^ d{hjUj) dZn 


(72) 


dx 


(74) 


5 The Navier-Stokes system 

Instead of considering the Euler system, we can also depart from the Navier-Stokes 
equations to derive a layer-averaged model. 

The model derivation is similar to what has been done in Section [3] for the Euler 
system. 


5.1 Layer averaging of the viscous terms 

In this paragraph and the both following, the components of the Cauchy stress tensor 
S are not specihed. It remains to hnd a layer-averaged formulation for the r.h.s. of 
Eq. (P, i.e. 


14 = 


+ 1/2 f dT^rr dE, 


^a-1/2 


+ 


+ 


d 


2 rr 


dx dz dx"^ 


Yj^xdzi 


as 


dx 


ZZ \ 1 

dz. 


We have 


dx 


Q r^a + l/2 


2a-1/2 


d n 

Sxi T I ^zxdzi S^z ) dz 


TSjjz |a-|-l/2 


dz, 


'a+l/2 


1-1/2 + 


dx 

dZa-\l2 


dx 


d n 

^xx T 7/ / ^zxd'Zi Sz 
dx J, 

d n 

Sxx T 7/ / Szxd^^i Sz 

dx L 


1 + 1/2 


1 - 1/2 
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In the expression Va we have the term 


g rz^+i /2 / Q rv \ 

— ^zxdzi dz 

\dxj, J 


ox \ ox J ^ Ox J ^ 

\ ^a-1/2 ^ ''^a + 1/2 


Y^yrdz 


dz^_ii2 n 


Tjzxdz 


' ^a-1/2 


d_ ^ p+i/2 ^ p 

/ Z^zxdz -\- 2^a+l/2 / ^zxdz 

OX \ OX I ^ , Ox J ^ 

\ ■'^a-l/2 ''^a + 1/2 

Za—1/2~Q f ^zxdz j , 

^a-l/2 / 


<92:„+i/2 ( d /■% 

O In/ 

ax \ax 


_ <9^a+l/2 <9 r ( dz^+i/2 

dx dxj, , V <92^ 

:< + l/2 -^^0 + 1/2 \ 


dz^_l/2 ( d ^ \ _dZa-ll2 dp ^ , fdZc,-l/2\ ^ 

dx Jz , dx dx A ^ ) S-U- 1 / 2 - 

5.2 Definitions and closure relation 

The expression of the viscous terms generally involving second order derivatives, their 
discretization requires quadrature formula that are not inherited from the layer-averaged 
discretization. In particular, at this step of the paper, we adopt the following notations 

^afe|„_,_i/2 ^ ^afe,0+1/2 ) 


and the following definitions. 


^ab\a ^ ^ab,a ) 


Yjabdz ~ ha^ab,a : 


with (a, b) G {x, z)"^. For the terms having the form 


'‘^0+1/2 


Z^abdz ^ 


a closure relation is needed and we choose the approximation 


'‘^a + l/2 10 — ZZ, 1 /9 

Ay ~ Q+1/2 0-1/2 y, _ U ^ Y 

~ o ^ab,a '^a^a-^ab^a' 

'a-l/2 
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For each interface we introduce the unit normal vector nQ_|_i /2 and the unit 

tangent vector tQ_|_i /2 given by: 


n 


a+l/2 


dz. 


+ 1/2 


dx 


Ca+1/2 


i + (^)' 

Then, for 0 ^ a ^ we have the following expression 


ta+l/2 


_ / Co+i/2 
’a+1/2 


to+1/2 ■ SQ,+i/2n„+i/2 — 


i + (%^) 

5^a+l/2 


\ ^CCX.Q-t-l/Z I 

OX \ ax 

which can be rewritten as 

ta+l/2 ■ Fla+i/2nQ,+i/2 = C^+;^/2'^a+l/2 ! 
by introducing the following notation, 


^xz^a-\-l/2 


+ 1/2 H-rj"*" ^ S^x,«+l/2 ~ S^2,a+l/2 ) ) ; (79) 


_ y. ^^0+1/2 , 

^q;+1/2 ^xz^ol-\-\I2 1 ^ 


{^xx,a+l/2 + 


dz. 


a+1/2, 


dx 


■^2:x,ct+l/2 


-^zz,a+l/2 


(80) 


(81) 


Remark that, for 0 ^ a ^ iV, the quantity tQ,+i /2 ■ FlQ,+i/ 2 nQ,+i /2 represents the tangential 
component of the stress tensors at the interface Zo,+i/ 2 - And for a = {0, A^}, the quan¬ 
tities (79) coincide with the boundary conditions and hence are given. More precisely 


(since Ci /2 = c?,) the Navier friction at bottom gives 




tl/2 ■ Fli/2ni/2 — —Ml — C’'1/2 Ci/2- 


(82) 


Compared to equation (10), velocity in the hrst layer ui is used since Ub is not a variable 


of our system. It is consistent with the convention (87) and dehnition (48). At the 
surface we have 

t7V+l/2 ■ '^N+l/2'^N+l/2 = ^N+1/2 (?n+1/2 = O' 


Remark 5.1 In (82) as in section^, we use the expression to consider a Navier 

friction at the bottom since on an impermeable boundary (10) is equivalent to For 
1 < a < A^ — 1, the flow can move across the interface +)a+i /2 o,nd we cannot give a 
formulation directly comparable to 


5.3 Layer-averaged Navier-Stokes system 

We have the following proposition. 
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Proposition 5.2 Using formula (77),(78) and (81), the layer-averaging applied to the 
Navier-Stokes system (|Tg-(|g completed with the boundary conditions (@-(g leads to 
the system 


N r. N 

o , o , 


dt ^ dx 
i=i j=i 

d . a 


(83) 


O . O f ^ 2 9 7 7 dzb „ „ 

-^{haUa) + ^ yhaU^ + -hallJ — ~9ha^^ + Wq,+i/2(jq;+1/2 ~ a-\l2 


A 

9a: 


^ U ( 

“l“7^ 1 ha^xx.a ha^zz,a “1“ ( haZa^zx,a 

'J ^a—1/2 


N 


+^a+l/2 Q^2 hj'Lzxj 

j=a+l 

+(^a+l/2 — (^Oi-1/2 


92 


N 




j=a 


1 d{haUa) ^ 9(/ljMj) ^ 9 Zq, 
UJa ~ ^ ~1~ Hq: , 


2 9x 


i=i 


9a: 


dx 


a = 1,..., N 


(84) 

(85) 


with the exchange terms Ga±i /2 given by (59) and the interface terms (Ta±i /2 given by 

(0. 

For smooth solutions, we obtain the balance 


m 


N 


9 


N 


Q:=l 

dZg 

dx 


<a=l 


9: 


^ ^ Ea I “1“ j ^ ^ Vjg I Eg F ^hgH hg{^xx,a ^zz,a) 


9 .1 


N 


haFjzx,a F hg Q^(^^hgYjzx,a F ^ ^ j '^f^Wghg'Fjzx, 


j=a+l 


N 


Q = 1 


N 


^ ^ ( dx ^ 22 , 0 ) 

a=l ^ 

'dWg dZgdUg 


F 


haF7^x,a F CTq,+i/2(Wo+1 — Ug 


dx dx dx 

yj^th Eg = E^F ^ ^ 

In (86), we use the convention 

Uq = Ui, Mat+i = Mat. 


n\ « 2 
V ) „3“l’ 

/ G 


( 86 ) 


(87) 


Before to give the proof of prop. 5.2, we make few comments concerning the layer¬ 


averaging of the Cauchy stress tensor components. 


Remark 5.3 Since the expression of the components of the Cauchy stress tensor are 


not specified, we are not able to precise all the terms in Eg. (86) and we only intend to 
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demonstrate that the energy balance (86) is consistent with (14). The nonnegativity of 


the right hand side of (86) has then to be verified when specifying the rheological model 


(as it is done below in the Newtonian case). 


Remark 5.4 After injecting the definition (81) of aa+ 1/2 in (86), it appears that the 
following terms in the right hand side of (86) 


N 




ol=1 


dx 


t) {fi^xx,a+l /2 ^22,0+1/2) ("Wa+l 


account for a layer-averaging of 


'Zb 


^ du 

zz'}dz, 


appearing in the right hand side of (14). Likewise, the term 




dx 


( 88 ) 


in the right hand side of (14) is discretized by 

I ^ \ , r ^ , dzadua 

^ \ ^xz,a+l/2[Ua+l Ua) + naZ^zx,a\ 


OL=l 


dz, 


a+1/2 

dx 


■^zx^ol-\-1/2 


(89) 

in the layer-average context of Eg. (86). A similar comparison can be done for the viscous 


terms involved in the left hand side of the two energy balances (14) and (86). 


Proof of proposition 5.2 The derivation of Egs 


been done to obtain the layer-averaged Euler system (72)-(74). Only the treatment of 
the viscous terms Va has to be specified. 


and (85) is similar to what has 


Using the definitions (77),(78), (81), for a = {l,iV} using the mimic of the boundary 
conditions it comes 


d ( d 

ka ~ I h(fTjxx a hojEzz,a 4“ 7^ / zTizxdz 

dx \ dx J ^ , 

\ •' Za_\/2 


d2 ^ Q2 

+Za+l/2-Q^ ^ hj'^zxj — Za-l/21^''^^hj^zx,j 
j=a+l 

+0'a+l/2 ~ Cra-1/2- 


N 


j=a 


The approximation (78) gives 


Va ~ Ra + C’‘«-|-l/2 ~ <^0-1/2 
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d 


2 2 

_ “ I V- , ^ I ^ /^a+1/2 “ ^a-1/2^ 

o I ^a^xx.OL ^a^zz.a ~r 1 ~ 

OX 


‘'a‘^xx,a "'a^zz,a ' \ 

OX V 

N 




+^a+l/2^^2 

j=a+l 
+<^0+1/2 ~ <Oa-l/2- 




N 


9^2 




J=Q 


For the energy balance we write 

-^Q^CK ( nrvh/y { ^ 

Ox 


a'^ay^xx^a ^zz,a J i V 

N 

dx 


, 5 /^^a+l/2 ^«-l/2^ 

+ MaT^ I -^^ 


2 

N 


zx^a. 


“I“'^a+l/2^a o ^ ^ hjYj;^x,j ^a—l/2^ao ^^hj^zx,^ 


a I ^xx,a ‘-‘zz 


-/l„( S 

A 

dx 


j=a+l 

s 




j=a 

2 2 
\dUa dUa d /^a+1/2 ~ ^0-1/2 , 


(9x 9x 

N 


{- 


N 


d . ^ 5 5 

(2:a+l/2'Wa)— E] + ^(^a-l/2'*^a)^ E/ 


dx 


dx 


j=a+l ]=a 

Notice that, using an integration by part, it comes that the three terms 

N 


a ( ^1+1/2- ^1-1/2 

9a; I “ 9a; V 2 


^zx,a'j N Z,x+l/2^a ^ ^ hj'^zx,j 

j=a+l 


9 


N 




V2'“a a E 


j=a 


appearing in Eg. (90) are a discretization of the guantity 


E 

dx 


Ur 


:«+i/2 Q n 

7^ / E^xdzidz J, 


^a-l/2 


dx 


in the energy balance Eg. (86). 
We can see that 


E 

dx 


9 


N 


Ua\Za+l/2Q^ E] 

j=a+l 


9 


AT 


,i ^01-1/2 


zx,j 


j=a 


9 


9 ^ 9 ^ 9 

(/Iq, F '^cf—1/2)7^ ^ ^ hj^zx,j Zqi—\I2~q^ ^ ^ hj^zx,j '^a—1/2 {hct^zx,a) 


j=a+l 


j=a+l 


dx 


N 


d f ( X ^ X y ^ (h T \ 

Un \ / j ^j^zx,j '^ 0 — 1 / 27 ^ \iT'aZ‘zx,aj 


j=a+l 
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Denoting RaUa the last three terms in Eg. (90), we write 


Ra^a I 

OX \ 


d A^a+1/2 ^a-l/2dUc 


dx 


^0+1/2 

~r ~ z Tir ^ 7 


N 


dx'^ 

N 


Q Q _^ Q Q _^ 

— T^(^n-+l/2WQ)^ hjTizxJ + -^{Za-l/2Ua)-^ hj'^zx,j 

j=a+l j=a 

2 ^2 


dx 

d_ 

dx 

( 

dx 


_ ^ f fU U U \\- ^ I ^«+l/2 ^a-ll2d‘^U^^ 

O I (^Q^a h(yk(yj^;zx,a 1 “1“ r\ ^ „ S 


dx'^ 


d d ^ d d 

— — (Zn+l/2Ua)-^ hjTizxJ + -^{Za-l/2Ua) 

j=a+l j=a 


N 


where (67) has been used. And simple manipulations give 

' dwa dZa dUa 


d 

RnUry = -— ( WnhryTl 


dx 

( 

dx 


Ja'^a^zx.cx 


d 


r\ ' o o I ha^zx^a r., {ha^zx,a) 

dx dx dx / dx 


N 


N 


d d 7 ^ d d X 

— — (Za+l/2Ua)-^ hjEj^xJ + ^ (^a-l/2'»a) ^ hjTizxJ 

j=a+l j=a 


d 1 

w cihcfTj 

dx ' 

+UJc 

d 

i+l/2^ 

dx 


dWa ^ dZa dUa 


hry^Z 


r\ I O O / ' ZX .CX. 

dx dx dx / 

N N 


d 


j=a+l 


,j Wa-ll2 hj^zx,ji 


j=a 


with Wa+i /2 defined by 


5(Zq,_|_i/2'Uq,) 5(Zq_|_i/2Mq,_|_i) 

^a+1/2 ka ^Q+1 


dx 


dx 


The two last terms of RaUa give a telescoping series and vanish when summing since 
Wi /2 = 0 and ^j^zx,j vanish when a = N. 

Finally, the guantity 


N 


'y ^ ^aV>ai 


Q=1 


gives the expression involving of the terms related to the Cauchy stress tensor in (86) 
proving the result. ■ 
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5.4 Newtonian fluids 

When considering a Newtonian fluid, the chosen form of the viscosity tensor is 


du 

dx' 

rdU dw y 

(93) 

dw 

, du dw y 

(94) 


where /r is a dynamic viscosity coefficient. 

When considering the fluid rheology is given by (93)-(94), thus leading to = —T^xx 
and Tixz = ^zx, prop. |5.2| becomes: 

Lemma 5.5 The layer-averaging applied to the Navier-Stokes system for a newtonian 
fluid gives 


r. N r. N 

o , o , 


dt 


(95) 


i=i i=i 

a . d 


o' ,, , o' /, 2 S' ; j dZb „ „ 

-^[haUa) + ^ \ haU^ + -ha-tl J — Ua+l/2(^a+l/2 ~ Ua-l/2(^a-l/2 


d 


dx 


d 


1 ‘^hoiTjxx.a ( haZofTiz^x.a 


dx 


32 ^ q2 

+Za+i/2-Q^ ^ ^ hjSj 

j=a+l 

+Cro+l/2 ~ (Za-l/2 1 

1 d{haUa) d{hjUj] 

- } -0 = 1,...,^ 


N 


Wry = -- 


2 dx 


i=i 


dx 


j=a 


dZr, 


dx ’ 


(96) 

(97) 


where exchange terms Ga±i /2 ooe still given by (59) and the interface terms (Ta±i /2 
defined by (81) are here reduced to 


_ OV- 5^a+l/2 , ^ 

^q+1/2 ^^xx,afl-l/2 “T ^zx^ol-\-1/2 

For smooth solutions, we obtain the balance 


1 - 


dZr 




dx / 


(98) 


m 


N 


d 


N 


^ ^ Ea I + j ^ ^ Uct I Ea + —haH — 2haTi 


Q = 1 

dZry 

dx 


ya=l 


^-1 

2 


d rl 


N 


ha'T‘zx,a F hry Q^i^^hrfFjzx^a F ^ ^ hj'Fj^x,j')j j ^ ^ WghryYj 


j=Oi+l 


N 


a=l 


N 

E 

a=l 


dUry^, _ . /'dWry dZadUa 


2h(yTjxx.a F 


+ 


^Oi XX.CX VO O O 

dx V ox dx dx 


jhryT^x,a F CrQ-|-i/2 (^«+l ^a) j 3^1’ 


(99) 
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If we look at the energy balance for the continuous setting (14), we have, by using 

m right hand side, 

+ ( 100 ) 

/X ^ ' Cz, 


(93)-(94), the following non-positive right hand side, 

n I 


of the layer-averaged model leads to 


N 


whereas, after including (98) in (99), the right hand side of the discrete energy balance 

dZa+l/2 , (dWo, , dZadu^ 

7!' 

( 101 ) 


a=l 


.du. 


Re ^ ^ I 2 ‘^^xx,a+l/2i,'^a+l 


Ur 


dx 


+ 


O ' O ^-]hcX‘zx,a 
OX ox Ox / 


+ 


+ S,,,„+i/2(n„+l - Ma) (^1 - • 


The aim of the next proposition is to mimic (105). 


Proposition 5.6 The layer-averaging, given in lemma |5.5[ is applied to the Navier- 
Stokes system for a newtonian fluid with the following consistent expressions of the 
rheology terms at the interface a -|- 1/2, 

^a-|-l/2^a;a;,Q:-|-l/2 hn-\-\^2'R‘zz,a+l/2 

= 21* f + 1 (.102) 


dx 


dx 


dx 


^(1+1/2“^ 2X,q:+1/2 


^a-\-l/2^xz,a-\-l/2 
= /i 


-fh U ( dWg+i dZa+ldUg+l ^ 

2\'^dx dx dx dx dx dx 


T('nQ,-|_i Ugf ( 1 


/ 92(q,+i/2\ 2\ 

\ dx ) ) 


(103) 


and, since the rheology terms are more related to elliptic than hyperbolic type behaviour, 
we used the centred approximation for the rheology terms at the layers a, 


Tjnh,g 


Safe,a+l/2 + 'R‘ab,a-l/2 


(104) 


with {a,b) G {x,zf. Then we obtain an energy ineguality since the right hand side of 


the discrete energy balance of the layer-averaged model, defined by (101), leads here to 

N 


Re = -Y^ 


hg+l/2 I ^2 , v^2 


a=0 


T 


xx,a+l/2 4” 0 + 1/2 j ^3^1 ' 


K 


(105) 
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Proof The expression (105) clearly mimics the continuous one given by (100). Moreover 
it is possible to exhibit a kind of consistency of the definitions (105)-(102). Indeed if we 
express the derivatives of the newtonian stress terms along the interface 0 + 1/2, on one 
hand, we have 


'du{x, Za+ll2{x, t),t) dZc,+l/2{x, t) 


= 2/i 


dx 


dx 


dzU{^X, Z,t)\z=z^^i/2{x,t'^ , 


which is consistent with (102). And, on the other hand, we have, 

^2a;|2=2^_,_i/2(a:,t) ~ T |2=2c«+l/2+,0 + ^x;W{x, Z,t)\z=z,y^i/2ix,t)^ ■ 

Additionally, we can write 

^ _dw{x,Zo,+xi2{A,f),t) 9z„+i/2(a;,t) 

dz,w{x,Z,t)\z=z^^^l 2 {x,t) dz'^yX, Z,t)\z=z„,+^/ 2 {x,t) 5 

and, using the incompressibility condition, we get, 

dzW{x, Z,t)\z=z,„j^l/2{x,t) ^x'^iAi ^A')\z=z„^j^i/2{x,t) ■ 

Therefore we have. 


O . dw{x,Za+l/ 2 ix,t)A) , 

dxW{x,Z,t)\z=z,„+,/2ix,t) = - ^ -+ 

dZa+l/2{x, t) f du{x, Za+l/2{x, t),t) dZc,+l/2{x, t) 


dx 


dx 


^^^)\z=Zc,+i/2ix,t)J 


Finally, this leads to the following expression 

'^zx\z=z,„+l/2(x,t) = T 


dw{x,Za+il2{x, t),t) ^ dZa+i/2{x,t) du{x, Za+x/2{x ,t) ,t) 


dx 


dx 

dzo,+i/2{x,t)- 

dx 


dx 


dzU{x, Z, t)|2=2^+i/2+,i) ) ) 


which is consistent with (103). 


The energy inequality is obtain by injecting (102), (103) and (104) in (101). 


Remark 5.7 We can remark in the lemma (5.5) that the rheology terms are both at the 
interface and in the layers. Thus an other strategy could be to defined them at the layer, 
and to average the terms at the interface. In this case, we have 


haA-‘xx,a h(fTi zz,a 


= 2fr [ h, 


dUa 
‘ dx 


^'^a+l/2 '^a+l ^ 1/2 '^a ^a—1 


dx 


dx 


(106) 
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a^zx,a 


hct^xz,a 

, dWa , , dZa dUa , Ma+1 “ “a , , 

I ^a~ - 1 “ ha~ -- 1 -:- I 1 — 


dx 


^CK —1 


dx dx 


dzr 


+1/2 Y 


dx / 


1 - 


dzr,- 


a-l/2 

dx 


( 107 ) 


which are also consistent expressions of the tensor, and the following averaging is intro¬ 
duced, 


■^ab^OL+ 1/2 


^a6,a+l “1“ ^a6,c 


(108) 


and leads to an energy ineguality, since the right hand side of the discrete energy balance 


of the layer-averaged model, defined by (101), leads here to 

N 

2 \ ^2 


a=l ' 


(109) 


This strategy seems to be more natural since, in the spirit of the layer-averaged model, the 
unknowns are mainly localised in the layers. However the main drawback is the stencil 
of the interface rheology terms which are not compact. For instance, the term 
will be expressed in function of Ua+ 2 ,Ua+i and u^-i. 


5.5 An extended Saint-Venant system 


In the simplified case of a single layer, the model given in prop. |5.2| corresponds to the 
classical Saint-Venant system but completed with rheology terms. 


Proposition 5.8 The classical Saint-Venant corresponds to the single-layer version of 
the layer-averaged Navier-Stokes system. With obvious notations, it is given by 

dH d 

w + 


¥^+(”’*l»’)-.»g 


d 


+ ( HTjxx ~ HTizz + 7 ^ ( 

dx \ dx \ 

_ ld{Hu) _d fH + 2zb 




dx^ ci 


For smooth solutions, we obtain the balance 


dE d /_/„ ^^,2 ^ \ 9 fd{H + 2zb)^ H d 


— + — M _ H(J:xx - ^zz) - 7 ^ 

dt dx\ \ 2 ^ 


dx V 


dx 


'^xz + ~J^{HTlxz) 

2 dx 


27 



















~ '^zz) 


„{dw I d{H + 2zi,) du\— k 
. dx^2 dx dx)^"" 


, 3 “ ’ 


wtthE = ^ + 1 (H + Zb) 


In the particular case of a newtonian fluid, the Saint-Venant system given in prop. 5^ 
reduces to 


2!|a.- -„g 

d [ du d f{H + ZbY-zldw 

— i — 

_ ld{Hu) _d fH + 2zb 
^ “ “2 dx V 2 


( 110 ) 




( 112 ) 


For smooth solutions, we obtain the energy balance 


dE d ( g 2 TT^'^\ ^ f fd{EI + 2zb) dw H d 


/i 


dt ClX ' ' f)rr.) 


dx^ dx 


dx 


dx ^ 2 dx^^dx 


i;r 7 T 2 


H dw 
2 dx 




du^ 


^,'dw ld{H + 2zb)duY k _2 

+ 2 ax & - 5 “' 


Remark 5.9 Notice that, compared to the classical viscous Saint-Venant system m, 


the model (110)-(113) has complementary terms. 


6 Conclusion 

We have proposed a layer-averaged discretization for the approximation of the incom¬ 
pressible free surface Euler and Navier-Stokes equations. The obtained models do not 
rely on any asymptotic expansion but on a criterion of minimal kinetic energy. Notice 
also that the layer averaging for the Navier-Stokes system has been carried out for a 
fluid with a general rheology. 

Since these models are formulated over a fixed domain, it is possible to derive efficient 
numerical techniques for their approximation. For the approximation of the proposed 
models, a finite volume strategy - relying on a kinetic interpretation and satisfying 
stability properties such as a fully discrete entropy inequality - will be published in a 
forthcoming paper. 
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